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Abstract: A variable annuity contract with Guaranteed Minimum Withdrawal Benefit (GMWB) promises to 
return the entire initial investment through cash withdrawals during the policy life plus the remaining account 
balance at maturity, regardless of the portfolio performance. We assume that market is complete in financial 
risk and also there is no mortality risk (in the event of policyholder death, the contract is maintained by bene¬ 
ficiary), thus the annuity price can be expressed as an appropriate expectation. Under the optimal withdrawal 
strategy of a policyholder, the pricing of variable annuities with GMWB is an optimal stochastic control prob¬ 
lem. The surrender feature available in marketed products allows termination of the contract before maturity, 
making it also an optimal stopping problem. 


Although the surrender feature is quite common in variable annuity contracts, there appears to be no published 
analysis and results for this feature in GMWB under optimal policyh older behavior - results found in the litera¬ 
ture so far are consistent with the absence of such a feature. Recently, Azimzadeh and Forsvth ( 20141) prove the 
existence of an optimal bang-bang control for a Guaranteed Lifelong Withdrawal Benefits (GLWB) contract. 
In particular, they find that the holder of a GLWB can maximize a writers losses by only ever performing non¬ 
withdrawal, withdrawal at exactly the contract rate, or full surrender. This dramatically reduces the optimal 
strategy space. However, they also demonstrate that the related GMWB contract is not convexity preserving. 



actually performed calculations with surrender option on top of optimal withdrawal strategy. Also, it is of 
practical interest to see how the much simpler bang-bang strategy, although not optimal for GMWB, compares 
with optimal GMWB strategy with surrender option. 

Recently, in lLuo and Shevchenkol(l2015al) . we have developed a new efficient numerical algorithm for pricing 
GMWB contracts in the case when transition density of the underlying asset between withdrawal dates or its 
moments are known. This algorithm relies on computing the expected contract value through a high order 
Gauss-Hermite quadrature applied on a cubic spline interpolation and much faster than the standard partial 
differential equation methods. In this paper we extend our algorithm to include surrender option in GMWB 
and compare prices under different policyholder strategies: optimal, static and bang-bang. Results indicate that 
following a simple but sub-optimal bang-bang strategy does not lead to significant reduction in the price or 
equivalently in the fee, in comparison with the optimal strategy. We also observed that the extra value added 
by the surrender option strongly depends on volatility and the penalty charge, among other factors such as 
contractual rate, maturity and interest rate etc. At high volatility or at low penalty charge, the surrender feature 
adds very significant value to the GMWB contract - the required fair fee is more than doubled in some cases; 
thus it is critical to account for surrender feature in pricing of real products. We also performed calculations for 
static withdrawal with surrender option, which is the same as bang-bang minus the “no-withdrawal” choice. 
We find that the fee for such contract is only less than 1% smaller when compared to the case of bang-bang 
strategy, meaning that the “no-withdrawal” option adds little value to the contract. 


Keywords: Variable annuity, optimal stochastic control, optimal stopping time, bang-bang control. Guaran¬ 
teed Minimum Withdrawal Benefit, surrender option, Gauss-Hermite quadrature, cubic spline. 
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1 INTRODUCTION 


The world population is becoming older fast. As a consequence the age-related spending is projected to rise 
dramatically in the coming decades in all the developed countries. Increasingly governments in the developed 
world realize they cannot afford paying sufficient public pensions and are looking for innovations in retirement 
income product market. In this paper we consider a variable annuity contract with Guaranteed Minimum 
Withdrawal Benefit (GMWB) with option to surrender the contract before maturity. This contract promises to 
return the entire initial investment through cash withdrawals during the policy life plus the remaining account 
balance at maturity, regardless of the portfolio performance. Thus even when the account of the policyholder 
falls to zero before maturity, GMWB feature will continue to provide the guaranteed cashflows. In addition, 
we allow the option to surrender the contract before the maturity which is a standard feature of real products on 
the market. GMWB allows the policyholder to withdraw funds below or at contractual rate without penalty and 
above the contractual rate with some penalty. If the policyholder behaves passively and the withdraw amount 
at each withdrawal date is predetermined at the beginning of the contract, then the behavior of the policyholder 
is called “static”. In this case the paths of the account can be simulated and a standard Monte Carlo simulation 
method can be used to price the GMWB. On the other hand if the policyholder optimally decide the amount 
of withdraw at each withdrawal date, then the behavior of the policyholder is called “dynamic”. Under the 
optimal withdrawal strategy of a policyholder, the pricing of variable annuities with GMWB becomes an 
optimal stochastic control problem; and adding surrender feature makes it also an optimal stopping problem. 


The variable annuities with GMWB feature under dyna mic and static withdrawal strateg i es have been consid 


ered in a numbe r of papers over the l ast de c ade, e.g. iMilevskv and Salisbury! (l200i 
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Dai et al. ( 2008ll. Huang and Forsyth ( 2012 ): Huang and KwokI ( 20141) . Bacinello et al 


Bauer et 51 ( 2008 ). 


(1201 Ih . Recently, 

Azimzadeh and ForsvthI (12014ft prove the existence of an optimal bang-bang control for a Guaranteed Life¬ 


long Withdrawal Benefits (GLWB) contract. In particular, they find that the holder of a GLWB can maximize 
a contract writer loss by only ever performing non-withdrawal, withdrawal at exactly the contract rate, or full 
surrender. This dramatically reduces the optimal strategy space. However, they also demonstrate that the re¬ 
lated GMWB contract is not convexity preserving, and hence does not satisfy the bang-bang principle other 
than in certai n degenerate case s. F or GMWB under optimal withdrawal assumption, the numerical algorithms 
developed by Dai et alJ ( 20081) and Chen and Forsvth (2008) appear to be the only ones found in the literature, 
and both are based on solving corresponding partial differential equation (PDF) via finite difference method. 


In the case when transition density of the underlying wealth process between withdrawal dates or its moments 
are known in closed form, often it can be more convenient and more efficient to utilize direct integration 
methods to calculate the re quired annuity value expectations i n backward time-stepping procedure. Such an 
algorithm was developed in Luo and Shevchenko ( 2014 , 2015 ah for solving optimal stochastic control problem 
in pricing GMWB variable annuity. This allows to get virtually instant results for typical GMWB annuity 
prices on the standard desktop PC. In this paper we adopt this algorithm to price variable annuities with 
GMWB with surrender option under static, dynamic, and simplified bang-bang withdrawal strategies. To our 
best knowledge, there are no publications presenting results for GMWB with both optimal withdrawal and 
surrender features. 


In the next section we describe the GMWB product with discrete withdrawals, the underlying stochastic model 
and the optimization problem. Section [3] describes numerical algorithm utilized for pricing In Section H] 
numerical results for the fair fees under a series GMWB contract conditions are presented. Concluding remarks 
are given in Section|5] 


2 Model 


We assume that market is complete in financial risk and also there is no mortality risk (in the event of policy¬ 
holder death, the contract is maintained by beneficiary), thus the annuity price can be expressed as expectation 
under the risk neutral process of underlying asset. Let S{t) denote the value of the reference portfolio of assets 
(mutual fund index, etc.) underlying the variable annuity policy at time t that under no-arbitrage condition 
follows the risk neutral stochastic process 

dS{t) = r{t)S{t)dt a{t)S{t)dB{t), (1) 

where B{t) is the standard Wiener process, r{t) is risk free interest rate and a{t) is volatility. For simplicity 
hereafter we assume that model parameters are piece-wise constant functions of time for time discretization 
0 = fo < G < ■ ■ ■ < tN = T, where to = 0 is today and T is annuity contract maturity. Denote correspond- 
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ing asset values as S{tQ),..., S{tN)', and risk-free interest rate and volatility as ri,... ,rN and ai,..., ctat 
respectively. That is, ri is the interest rate for time teriod {to, <i]; r 2 is for (fi; t 2 ], etc and similar for volatility. 

The premium paid by policyholder upfront at to is invested into the reference portfolio of risky assets S{t). 
Denote the value of this variable annuity account (hereafter referred to as wealth account) at time f as VL(f), 
i.e. the upfront premium paid by policyholder is W{0). GMWB guarantees the return of the premium via 
withdrawals 7 „ > 0 allowed at times n = 1,... ,N. Let N.^ denote the number of withdrawals in a 
year (e.g. N.^ = 12 for a monthly withdrawal), then the total number of withdrawals N = \ x T ]. 
The total of withdrawals cannot exceed the guarantee W (0) and withdrawals can be different from contractual 
(guaranteed) withdrawal Gn = hL(0)(f„ — tn-i)/T, with penalties imposed if 7 ^ > Gn. Denote the annual 
contractual rate as p = 1/T. 

Denote the value of the guarantee at time t as A{t), hereafter referred to as guarantee account. Obviously, 
A(0) = W{0). For clarity of notation, denote the time immediately before t (i.e. before withdrawal) as t~, 
and immediately after t (i.e. after withdrawal) as Then the guarantee balance evolves as 

A(f+)=A(f-)- 7 „ = A(f+_i)- 7 „, n=l,2,...,iV (2) 

with A(T+) = 0, i.e. W (0) = A(0) > 71 + • • • + 7 Ar and > J2k=n T*- ^he account balance A{t) 

remains unchanged within the interval {tn-i, tn), n = 1,2,..., N. 

In the case of reference portfolio process O, the wealth account W (t) evolves as 
W{t-) = 

W{t+) = max(VL(f;;;)-7„,0) , n = l,2,...,N, 


(3) 

(4) 


where = tn — tn-i, Zn are iid standard Normal random variables and a is the annual fee charged by the 
insurance company. If the account balance becomes zero or negative, then it will stay zero till maturity. 

The cashflow received by the policyholder at withdrawal time is given by 


r 7n. if 0 < 7„ < Gn, 

\ Gn + (1 ~ P)i'ln — Gn), if "fn > Gn, 


(5) 


where is contractual withdrawal. That is, penalty is applied if withdrawal 7 ^ exceeds G„, i.e. /3 S [0,1] is 
the penalty applied to the portion of withdrawal above G„. 

If the policyholder decides to surrender at time slice t G (1,..., — 1), then policyholder receives cashflow 

Dt{W{G), A{tr)) and contract stops. For numerical example we assume that 


Dr{W{tr),A{tr)) ■= Gr{uiax{W{tr),A{tr))); 


(6) 


other standard surrender conditions can easily be implemented. 

Denote the value of variable annuity at time t as Qt{W {t), A{t)), i.e. it is determined by values of the wealth 
and guarantee accounts W{t) and A{t). At maturity, if not surrendered earlier, the policyholder takes the 
maximum between the remaining guarantee withdrawal net of penalty charge and the remaining balance of the 
personal account, i.e. the final payoff is 

Qt-)W{T-),A{T-)) := hM{W{T-),A{T-)) = max{W{T-),GN{A{T-))) . (7) 


Under the above assumptions/conditions, the fair no-arbitrage value of the annuity at time to is 


Qto (W(to),A(to)) = max Et^ 


B(0,r)Dn(W(t-),A(t-))If,^<T} 


N 


+B{0, N)hN{W{T ),A{T ))(1 — -f E^(0 g-)Q(7,) 


N = mln{T,N) -1, (8) 


where Bq „ = exp(— Jj" r{T)dT) is discounting factor and Ipj. is indicator function. Note that the today’s 
value of the annuity policy (5o(W^(0), ^(0)) is a function of policy fee a. Here, r is stopping time and 
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7 i,..., 7 Ar_i are the control variables chosen to maximize the expected value of discounted cashflows, and 
expectation Eo[-] is taken under the risk-neutral process conditional on Wq and Aq. The fair fee value of a 
corresponds to Qo {W (0), ^(0)) = kL(0). It is important to note that control variables and stopping time can 
be different for different realizations of underlying process and moreover the control variable 7 „ affects the 
transition law of the underlying wealth process from to tn+i- Overall, evaluating GMWB with surrender 
feature is solving optimal stochastic control problem with optimal stopping. 

Denote the state vector at time as Xn = A{t~)). Given that X = (Xi,... ,Xpf) is Markov 

process, it is easy to recognize that the annuity valuation under the optimal withdrawal strategy (O is optimal 
stochastic control problem for Markov process that can be solved recursively to find annuity value Qt^ (x) at 
tn, n = N — 1,... ,0 via backward induction 


Qt„ {x) = max sup 


Cni'JniXn)) + 6 




Qtr,+i{x')Kt^{dx'\x,jn) ] ,Dn{x) (9) 


Starting from final condition Qt{x) = max (VL(T“), C'Ar(A(T“))). Here Kt^{dx'\x,"fn) is the stochas¬ 
tic kernel representing probability to reach state in dx' at time tn+i if the withdrawal (action) 7 „ is ap- 
plied in the state x at tim e For a good textbook treatment of stochastic control problem in finance, see 
Bauerle and Rieden (1201 ih . Explicitly, this backward recursion can be solved as follows. 

The annuity price at any time t for a fixed A{t) is a function of W only. Note A{ti^_^) = A{t~) = A 
is constant in the period f“). Thus in a backward time-stepping setting (similar to a finite difference 

scheme) the annuity value at time t = can be evaluated as the following expectation 


{W{t+_,),A)=Et., 


(Wit-), A) \Wit+_,),A 


( 10 ) 


Assuming the conditional probability distribution density of W{t^) given is known as 

Pn{witn)\witn-i)), then the above expectation can be evaluated by 

p+oo 

(lL(<+_i),A) = / e-^-‘^^-pniw\Witt_,))Q,-iw,A)dw. (11) 

Jo 

In the case of wealth process (O the transition density pn{witn)\witn-i)) is known in closed form and we 
will use Gauss-Hermite quadrature for the evaluation of the above integration over an infinite domain. The 
required continuous function QtiW, A) will be approximated by a cubic spline interpolation on a discretized 
grid in the W space. 

Any change of A(f) only occurs at withdrawal dates. After the amount 7 „ is drawn at the wealth account 
reduces from to kF(f+) = max(kF(f“) — 7 n, 0 ), and the guarantee balance drops from A(f“) to 

^(^n) = ~ In- Thus the jump Condition of (5t (IT, A) across is given by 


Qt;:iWit-),A{t-)) 

= max[ max [Q^+{inaxiW{t-) - ^ri,0), A{t-) - jn) + Cni^n)], Dr,{Wit-), A{t-))] . (12) 


For optimal strategy, we chose a value for 7 „ under the restriction 0 < 7 „ < A(f„ ) to maximize the function 
value Q^- [W, A) in (fT^ . Repeatedly applying (fTTll and (fl^ backwards in time starting from 

Qt-(W"(r-),A(T-)) =max(lL(r-),CAr(A(T-))) (13) 

gives us annuity value at f = 0. 

In additional to dynamic and static strategies, in this paper we also consider bang-bang strategy which is 
simplified sup-optimal strategy where the policy holder at each tn can either make withdrawal at contractual 
rate G„, make no withdrawal or surrender. 


3 Numerical algorithm 

A very detailed descr i ption o f the algorithm that we adapt for pricing GMWB with surrender can be found in 
Luo and Shevchenko ( 2015al) . Below we outline the main steps. We discretize the asset domain [Wmin) hEmax] 
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by Wmin = Wq < Wi ,..., Wm = Wmax , where Wmin and ILmax are the lower and upper boundary, 
respectively. The idea is to find annuity values at all these grid points at each time step from t~ to 
through integration (fTTTi . starting at maturity t = = T~. At each time step we evaluate the integration (fTTT) 

for every grid point by a high accuracy Gauss-Hermite numerical quadrature; it can also be accomplished by 
solving corresponding PDE using finite difference method that we implemented for benchmarking. At time 
step t~ t^-i, the annuity value at t = t~ is known only at grid points Wm, m = 0,1,..., M. In order 
to approximate the continuous function Qt{W, A) from the values at the discrete grid points, we use the cubic 
spline interpolation which is smooth in the first derivative and continuous in the second derivative. 

For guarantee account balance variable A, we introduce an auxiliary finite grid 0 = Ai < • • • < Aj = W (0) 
to track the remaining guarantee balance A, where J is the total number of nodes in the guarantee balance 
amount coordinate. For each Aj, we associate a continuous solution QtiW, Aj). At every jump we let A to 
be one of the grid points Aj, 1 < j < J. For any W = Wm, m = 0,1,..., M and A = Aj, j = 1,... ,J , 
given that withdrawal amount can only take the pre-defined values 'y = Aj — Ak, k = 1,2,j, irrespective 
of time tn and account value Wm, the jump condition (IT2l) takes the following form for the specific numerical 
setting 

Qt-{Wm,Aj) = max ( max [Q +(max(lLm - Aj + Ak,0),Ak) + Cn{Aj - Ak)], Dn{Wm, Aj)] . (14) 
" Vi<fc<i " / 

Overall we have J numerical solutions (obtained through integration) to track, corresponding to each of the 
Aj value, 1 < j < J. Stepping backward in time, we find Qq{W{Q), 4l(0)) that depends on the policy fee a. 
Finally, we calculate fair fee value of a corresponding to Qo(W^(0), ^(0)) = W{A) that obviously requires 
iterative process. 

4 Numerical Results 


Below we present numerical results for fair fee of GMWB with surrender option under optimal and suboptimal 
bang-bang withdrawal strategies. For convenience we denote results for optimal withdr awal strategy without 
surrend er option as GMWB, and with surrender option as GMWB-S. As discussed in Luo and Shevchenkol 
( 2015ah . only very few results for GMWB under dynamic policyholder behavior can be found in the literature, 
and these results are for G MWB without the surren der option. For validation purposes, perhaps the most 
accurate results are found in ichen and ForsvthI ( 20081) . which were obtained with a very fine mesh in a detailed 
convergence study . As s hown in Table [T] our GMWB results for fair fee compare very well with those of 


Chen and ForsvthI (l2008ll . The maximum absolute difference in the fair fee rate between the two numerical 


studies is only 0.3 basis point (a basis point is 0.01% of a rate). 


Table[T]shows some very interesting comparison among GMWB, GMWB-S and bang-bang results. At volatil¬ 
ity a = 0.2, the fair fee for GMWB-S is virtually the same as GMWB, meaning surrender adds little value to 
the optimal strategy; at high volatility cr = 0.3, fees for GMWB-S are significantly higher than GMWB, up 
to 50% higher at the half-yearly frequency. This may suggest that at high volatility it is optimal to surrender 
at high values of account balance or guarantee level. In addition, it also shows higher frequency gives higher 
extra value to the surrender option. Comparing bang bang with GMWB-S, the fees are below the optimal 
strategy as expected, but the values are not very significantly lower at both volatility values - it is only about 
10% lower at most. 


Figure[T]shows curves of the fee as a function of the contractual annual withdraw rate, given cr = 0.2, r = 0.05 
and /3 = 0.1. It compares four cases: static (without surrender), GMWB, GMWB-S and bang-bang, all at 
quarterly withdrawal frequency with 10% penalty charge, i.e. /3 = 0.1. This comparison also shows GMWB- 
S and GMWB have virtually the same fees at cr = 0.2, and bang-bnag is only slightly below GMWB-S, 
confirming results in Table[T] However, at the same volatility cr = 0.2, new features emerge if we reduce the 
penalty charge from /3 = 0.1 to ^ = 0.05, as shown in Figure |3 When the penalty charge is reduced and 
all other parameters are unchanged, the surrender option adds more significant value to GMWB - in fact the 
fees are more than doubled at low to moderate contractual withdraw rate (or equivalently long or moderate 
maturity), i.e. fees for GMWB-S are more than twice as those for GMWB. With reduced penalty, fees for 
bang-bang are still close to the optimal strategy with surrender option, the GMWB-S. 

We also performed calculations for static withdrawal with surrender option, which is the same as bang-bang 
minus the “no-withdrawal” choice. We find the fee for such contract is only less than 1% smaller than the 
bang-bang strategy, meaning the “no-withdrawal” option adds little value to the contract. Finally, different 
penalty functions can be applied to the surrender (i.e. surrender cashflow can be different from (|6ll). For 
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example, instead of penalizing only the amount exceeding the contractual withdrawal rate, we can penalize the 
entire termination amount. In this case we find both GMWB-S and bang-bang yield only slightly lower fees for 
a given /3 - this is perhaps not very surprising since when it is optimal to surrender, the amount must be much 
higher than the contractual rate, thus penalizing the entire amount is not much more severe than penalizing 
only the exceeded part. 


frequency 

volatility 

Chen & Forsyth 

GMWB 

GMWB-S 

Bang Bang 

yearly 

0.2 

129.1 

129.1 

129.2 

123.9 

half-yearly 

0.2 

133.5 

133.7 

134.0 

125.6 

yearly 

0.3 

293.3 

293.5 

418.4 

392.9 

half-yearly 

0.3 

302.4 

302.7 

456.5 

410.7 


Table 1. Comparison of fair fee a in basis points (a basis point is 0.01%) between results of GMWB, GMWB- 
S and bang-bang. Results under “Chen & Forsyth” are for GMWB. The input parameters are g = 10%, 
P = 10%, r = 5% and cr = 0.2. The withdrawal frequency is quarterly. 



g 


Figure 1. Fair fee a as a function of annual guarantee rate g for static, GMWB, GMWB-S and bang-bang at 
a quarterly withdraw rate. The fixed input parameters are /3 = 10%, r = 5% and cr = 0.2. 



Figure 2. Fair fee a as a function of annual guarantee rate g for static, GMWB, GMWB-S and bang-bang at 
a quarterly withdraw rate. The fixed input parameters are /3 = 5%, r = 5% and a = 0.2. 
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5 CONCLUSIONS 


In this paper we have developed numerical valuation of variable annuities with GMWB and surrender features 
under both static, dynamic (optimal) and bang-bang policyholder strategies. The results indicate that following 
a simple bang-bang strategy does not lead to significant reduction in the price or equivalently in the fee. We 
also observed that the extra value added by the surrender option very much depends on volatility and the 
penalty charge, among other facts such as contractual rate and maturity. At high volatility or at low penalty 
charge, the surrender feature adds very significant value to the GMWB contract - more than doubling in some 
cases; highlighting the importance of accounting for surrender feature in pricing of real products. We have 
assumed the policyholder will always live beyond the maturity date or there is always someone there to make 
optimal withdrawal decisions for the entire duration of the contract. It is not difficult to consider adding some 
death benefi ts on top of GMWB, i . e. comb ining GMWB with some kind of life insurance as it is done in our 
recent paper Luo and Shevchen^ ( 2015^ considering both market process and death process. Further work 
includes admitting other stochastic risk factors such as stochastic interest rate or volatility. 


Acknowledgement 


We gratefully acknowledge financial support by the CSIRO-Monash Superannuation Research Cluster, a col¬ 
laboration among CSIRO, Monash University, Griffith University, the University of Western Australia, the 

University of Warwick, and stakeholders of the retirement system in the interest of better outcomes for all. 

References 

Azimzadeh, Y. and P. A. Forsyth (2014). The existence of optimal bang-bang controls for gmxb contracts. 
Working paper of University of Waterloo. 

Bacinello, A., P. Millossovich, A. Olivieri, and E. Pitacco (2011). Variable annuities: a unifying valuation 
approach. Insurance: Mathematics and Economics 49(1), 285-297. 

Bauer, D., A. Kling, and J. Russ (2008). A universal pricing framework for guaranteed minimum benefits in 
variable annuities. ASTINBulletin 38(2), 621-651. 

Bauerle, N. and U. Rieder (2011). Markov Decision Processes with Applications to Finance. Springer, Berlin. 

Chen, Z. and P. Forsyth (2008). A numerical scheme for the impulse control formulation for pricing variable 
annuities with a guaranteed minimum withdrawal benefit (gmwb). Numerische Mathematik 109(A), 535- 
569. 

Dai, M., K. Y. Kwok, and J. Zong (2008). Guaranteed minimum withdrawal benefit in variable annuities. 
Mathematical Finance 18(A), 595-611. 

Huang, Y. and P. A. Forsyth (2012). Analysis of a penalty method for pricing a guaranteed minimum with¬ 
drawal benefit (GMWB). Journal of Numerical Analysis 32, 320-351. 

Huang, Y. and Y. K. Kwok (2014). Analysis of optimal dynamic withdrawal policies in withdrawal guarantee 
products. Journal of Economic Dynamics and Control 45, 19^3. 

Luo, X. and P. V. Shevchenko (2014). Fast and simple method for pricing exotic options using gauss- 
hermite quadrature on a cubic spline interpolation. Journal of Financial Engineering 1(A). DOT 
10.1142/S2345768614500330. 

Luo, X. and P. V. Shevchenko (2015a). Fast numerical method for pricing of variable annuities with guar¬ 
anteed minimum withdrawal benefit under optimal withdrawal strategy. To appear in Journal of Financial 
Engineering, preprint is available at http://arxiv.org/abs/1410.8609. 

Luo, X. and P. V. Shevchenko (2015b). Valuation of variable annuities with guaranteed minimum withdrawal 
and death benefits via stochastic control optimization. Insurance: Mathematics and Economics 62, 5-15. 

Milevsky, M. A. and T. S. Salisbury (2006). Financial valuation of guaranteed minimum withdrawal benefits. 
Insurance: Mathematics and Economics 38(1), 21-38. 








